Quantum gray solitons in confining potentials 
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We define and study hole-like excitations (the Lieb II mode) in a weakly interacting Bose liquid 
subject to external confinement. These excitations are obtained by semiclassical quantization of gray 
solitons propagating on top of a Thomas- Fermi background. Radiation of phonons by an accelerated 
gray soliton leads to a finite life-time of these excitations. It is shown that, for a large number of 
trapped atoms, most of the Lieb II levels can be experimentally resolved. 



Ultra cold atoms restricted to move in one dimension 
(Id) are now routinely used for experimental realization 
and investigation of the physics of strongly interacting 
many-body systems. Due to advances in measurement 
techniques, the dynamical response of one-dimensional 
quantum gases has recently gained a central role in these 
studies pJi3j. 

Dynamics is ultimately related to the nature of excita- 
tion spectrum. The latter is given by the Bethe Ansatz 
solution of Lieb-Liniger model [4 for Id ultra cold bosons 
with short-range interactions [5|i6j . As shown by Lieb [7] 
the excitations can be decomposed into a superposition of 
particle-like (Lieb I) and hole-like (Lieb II) excitations in 
one-to-one correspondence with elementary excitations of 
one-dimensional fermions. For weak coupling, the disper- 
sion of Lieb I excitations can be obtained semiclassically 
by linearization of the Gross-Pitaevskii Equation (GPE) 
resulting in phonon-like spectrum. Lieb I excitations de- 
termine the thermodynamic properties of the system and 
control the power-law decay of correlation functions. 

Later, it was shown by Kulish, Manakov and Faddeev 
[8] that the Lieb II mode can be associated with another 
semiclassical object - the gray soliton [9 . The latter is a 
hole-like solution of GPE representing localized density 
and supercurrent dip propagating on top a of constant 
background. It is a mean field picture of many quantum 
Lieb II excitations forming a coherent wavepacket. This 
classical description holds as long as quantum fluctua- 
tions remain small [10 due to a large negative effective 
mass resulting from the many particles expelled from the 
soliton core . The same large parameter is responsible for 
the exponentially suppressed probability to create gray 
solitons using local perturbations [11 . 

Even a shallow axial trapping potential present in the 
experiments breaks the integrability of the Lieb-Liniger 
model and the Lieb classification of the excitations is 
apparently lost. While Lieb I phonons can still be defined 
by linearization of GPE around the appropriate density 
profile and quantizing the resulting collective oscillations 
[m |13], no analogous procedure exists for the Lieb II 
mode. Classically, however, gray solitons do survive the 
harmonic confinement if the latter is sufficiently smooth 




FIG. 1. Schematic picture of trapped Lieb II modes obtained 
by quantizing motion of gray solitons in the non-uniform back- 
ground. The non- uniform density profile leads to the finite 
life-time r and precludes resolution of low-lying levels. 



on the length scale of the soliton [TH [15]. They have 
recently been observed in several experiments p^BHTS] in 
agreement with theoretical studies [I9l [20] . The physics 
of underlying quantum states has not been addressed yet. 

The aim of this Letter is to define and study the finite- 
size analogue of the Lieb II mode in a trapped system as 
corresponding quantum excitations. This is done by ap- 
plying Bohr-Sommerfeld quantization rule to the classical 
dynamics of gray soliton. The obtained quasiparticle lev- 
els, shown schematically in Fig. [l] have a finite life-time 
due to interactions of the soliton with the background ex- 
citations. The standard mechanism of soliton decay due 
to the scattering of thermal phonons [19 cannot be used 
here as it becomes inefficient for the exactly integrable 
GPE [21]. Even for non-integrable interactions the life- 
time diverges as inverse fourth power of temperature T 
[22] and solitons propagating in uniform background are 
absolutely stable at T = due to conservation of energy 
and momentum. The non-uniform density profile relaxes 
momentum conservation and leads to soliton decay even 
at T = 0. We calculate this fundamental limit on the 
life-time of the trapped quantum solitons. 
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Our main findings are as follows. For a system of A^tot 
interacting bosons trapped in a potential U (x) there are 
exactly Atot quantum states below the energy E^^ corre- 
sponding to the energy of stationary classical dark soliton 
localized at the maximum of the trapped density profile. 
We associate these states, shown schematically in Fig. [l] 
with the trapped Lieb II mode. In the case of harmonic 
confinement U{x) = we use the Thomas-Fermi 

density profile to show that £^ds = / V^)Ntot and the 
energies of the trapped Lieb II mode are given by the 
descending ladder 



En — En 



(1) 



where n = 0, 1, ... , Atot — 1- The energy of the high- 
est quantum state n = is reduced by zero-point os- 
cillations {hu/V2)/2 from the classical value £^ds- It is 
interesting to compare this result with the case of ex- 
tremely strong interactions, so-called Tonks-Girardeau 
limit [23], where hard-core bosons can be mapped into 
non-interacting fermions. In this case the Lieb II mode 
is obtained by creating holes in the filled Fermi sea of 
A^tot particles occupying energies En = E-p — huj{n-\- 1/2)^ 
below the Fermi energy Eip = ficjAtot- 

The trapped Lieb II states are not eigenstates but 
rather quasiparticle resonances with finite life- time. The 
life-time for a state with energy E is found to be 

8/i 



(2) 



where ja is the chemical potential of the condensate, 
= E/E^^, and the function F{x) is defined in Eq. (24). 
For high energies, E ~ £^ds7 we find a logarithmi- 
cally large life-time uor {Aji/hw) log [6£^ds/(£^ds — E)]^ 
while for low energy states, E <C £^ds we have lot 
{Sjii/3Jiw){E /E(^s)- The last expression defines an en- 
ergy scale = {fux)/ii)E^^ below which uor < 1 and 
Lieb II quasiparticle states cannot be resolved. This 
condition coincides with the classical picture of an over- 
damped soliton decaying before completing one period 
of oscillations in the trap. Another restriction, E ^ 
= E^^jK^I^ = /i/K^/^, arises from the requirement 
of the validity of the semiclassical treatment [11 . Here 
K = TT^ds//^ = {7rhuj/V2/jj)Ntot is the Luttinger param- 
eter in the center of the trap. For the experiment in 
Ref. [13 one has hw / ji ^ 10~^ and K ^ 10^, therefore 
£'* 0.01£^ds > E** and Eq. ([T]) describes accurately 
the energies of the most of the trapped Lieb II states 
and most of these are well defined quasiparticles, hence 
our semiclassical quantization of the soliton is justified a 
posteriori. 

Solitons in uniform background. — To derive the above 
results we consider the standard Lagrangian density for 
weakly interacting Id bosons C = ihi/jdt^p — Si^p^ip) with 
the energy density 
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2m 



(3) 



Here ilj{x^t)^ilj{x,t) are bosonic fields, m is the particle 
mass and ^ > characterizes repulsive interactions be- 
tween particles. In ([3| we have subtracted the constant 
contribution of the static background density n (fixed by 
the chemical potential ja = gn). Variation with respect 
to the fields leads to GPE, 



(4) 



As shown in [9^ it allows for a one-parameter family of 
gray solitons V^s(^ — ^(^)) well localized around the po- 
sition X moving with a constant velocity, X = Xq + Vt. 
According to [9 the form of these solutions is 



il)^{x)/^/n = cos 



^ . . ^ , f mgN 

z sm — tanh — z^^x 

2 2 \ 2h? 



(5) 



Parameters ^ and A" are integral characteristics of the 
soliton [24] representing the total phase drop across the 
soliton and the number of particles expelled from the 
soliton vicinity. The are both fixed by the velocity V of 
the soliton. 



^ = 2 arctan 



y2 



V 



N 



= Jdx(n- = y\/c2-y2. 



(6) 
(7) 



Here c = ^J gn/m is the sound velocity limiting the ve- 
locity of the soliton, V'^ < (? . 

Following Ref. [15' we wish to establish a description of 
a soliton as an effective particle. To this end we substi- 
tute the solitonic solution given by Eq. ([5| into Eq. (|3| 
and integrate it over the length of the system to obtain 
the energy of the soliton. 



4 hm 
3T 



2x3/2 _ mg 



2^3 



The canonical momentum of the soliton is 
Ps hn^ - mNV 



(8) 



(9) 



The contribution hn^ comes from a small background su- 
percurrent [25] which carries no energy but must be intro- 
duced to compensate for the phase drop ^. The second 
contribution to the momentum describes the deficit of 
A" particles moving with velocity V. Such an expression 
for the momentum is completely general and indepen- 
dent of the details of interactions between the particles 
[26] . In contrast, the energy E^ does depend on the form 
of interactions and the fact that it only depends on the 
particle deficit A" is a direct consequence of the quartic 
interaction of the CP energy density in Eq. (|3|. 

For our purposes it is convenient to use the canonical 
formalism where the momentum P, rather than velocity 
V serves as a control parameter. By formally inverting 
the relation ([9| we obtain the velocity as a function of 
momentum. 



X = V(P,n) 



(10) 
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Substituting V{P,n) into Eqs. (|6|, ^ and Eq. ^ yields 
$(P,n), N{P,n) and E^{P,n). The latter defines the 
dispersion of the soliton, or its effective Hamiltonian. In- 



deed, Eq. (10) constitutes the Hamiltonian equation of 
motion as we have V = {dE/dV) /{dP/dV) = dE/dP 
following from Eqs. (|8| and ([9|. Another Hamiltonian 
equation of motion states the conservation of momentum 
P = expected from the translational invariance. 

Solitons in Thomas- Fermi density profile. — 
Translational invariance is broken by the presence 
of external trapping potential U{x). In this case the 
equation of motion for the momentum will be altered. It 
should be stressed, however, that the external potential 
does not act on the coordinate X of the soliton directly, 
but rather on the particles of the Bose gas forming 
the soliton. To describe this situation and derive the 
equations of motion we follow the method pioneered 
in Ref.[T5]. We assume the non-uniform background 
density is given by the Thomas-Fermi density profile 



gn{x) = iJi — U{x). 



(11) 



For symmetric U{x), a nonzero density requires |x| < 
where the Thomas-Fermi radius R is found from U {R) = 
/i. The local sound velocity is c{x) = ^J gn{x)/m. 

For a sufficiently smooth potential and large R the 
density n{x) changes smoothly on the typical length scale 
of the soliton. In this case one can substitute n ti{X) 
into expressions for $(P, n), N{P^n) and E{P^n). The 
latter defines the effective soliton Hamiltonian in the non- 
uniform background. 



H{P,X)=E{P,n{X)), 



(12) 



which depends on the external potential U{X) only via 
the Thomas- Fermi density profile (11). 



Classical and quantum dynamics of gray solitons. — 
The Hamiltonian ([T2| generates the evolution for any 
dynamic observable 0(P, X) through the corresponding 
Poisson bracket O = {0,H} = dxOdpH - dpOdxH. 
In particular, the number of expelled particles N is con- 
served during the dynamics as follows from Eq. ([7|). The 
conserved combination with dimensions of energy 



mV'^ - mc^{X) = mX^ + U{x) - /i 



(13) 



can be used to solve for the dynamics of the soliton by 
mapping onto an effective particle with mass 2m moving 
in the potential U {X) as was done in Ref. ^ . Due to the 
positiveness of c'^{X) — V'^ the expression in Eq. (13) is 



always negative which restricts the motion to the region 
where U{X) < /i, i.e. inside the atomic cloud X < R. 

In what follows the potential U (x) is assumed to have 
only one minimum situated at X = 0. A classical point- 
like trajectory ^ = and X = describes the stationary 
dark soliton with ^ = ±7r at the maximum of the density 
profile corresponding to maximum energy £^ds- Any tra- 
jectory with non-zero velocity V has lower energy, thus 
the classical energy is bounded between and ^ds- 




FIG. 2. The classical phase space trajectories of a trapped 
gray soliton for energies between ^ = (horizontal line for 
P = 0) and E = Pds (points at X = and P = ±7rno). The 
shaded area defines the value of action variable /(P), Eq. ( p4| 
which is shown as function of energy to the right. The dashed 
line corresponds to an arbitrary trap potential U{x) while the 
straight solid line corresponds to an harmonic trap. Circles 
denote quantized values of 1(E) /h — n + 1/2 given by the 
Bohr-Sommerfeld rule. 



The typical phase-space portrait is shown in Fig. [2j 
As the velocity V changes sign at the turning points, the 
phase ^ given by Eq. (|6| undergoes a "umklapp" process 
and changes by 27r which in turn results in a correspond- 
ing jump of the momentum in Eq. ([9|. Calculating the 
oriented area enclosed by a given orbit in the phase space 
defines the action variable, 

I{E) = ^jp{E,X)dX. (14) 

as shown by shaded region in Fig. [2j Due to the hole- 
like character of gray soliton excitations, the action is a 
decreasing function of the energy. The maximum value of 



the action variable in Eq. ( 14 ) corresponds to the singular 



classical orbit which encloses the entire available phase 
space and on its two branches <I> = 0, 27r, V = ±c and 
N = 0. For this orbit the action variable takes the form 



In 



J dXhn{X) = HNu 



(15) 



The minimum value, / = is realized for the point-like 
trajectory with E = P^s- For intermediate energies the 



integral in Eq. ( 14 ) is a smooth function decreasing from 
/(O) = hNtot to /(Pds) = 0- The classical frequency of 
oscillations is given by dE/dl = Q{E) which is negative 
as expected for the hole-like excitations. 

Applying the standard Bohr-Sommerfeld semiclassical 
quantization rule, the energy levels are found from the 
condition I{En)/h = n + 1/2. Therefore there are exactly 
Xtot quantum states, for any trapping potential U{x). 
In the case of an harmonic confining potential U (X) = 
mcj^X^/2, the effective energy in Eq. ([l3| describes har- 
monic motion with energy-independent oscillation fre- 
quency ft = -cj/x/2 [HITS], hence Pds = {huo/ V2)Ntot 
and we obtain equidistant spectrum given by Eq. ([T]). 
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Interactions with phonons. — The picture presented 
above of solitons as Landau quasiparticles must be com- 
plemented with their interactions with the phonons. In- 
deed, a sohton propagating in a smooth non-uniform den- 
sity profile is analogous to an electron moving in static 
electric field created by some external sources. It is well 
known that its deceleration leads to radiation of electro- 
magnetic waves removing momentum and energy from 
the electron, known as Bremsstrahlung [27 . In the case of 
the gray soliton the emitted phonons lead to an eventual 
loss of energy and momentum from the soliton resulting 
in its acceleration due to its negative effective mass. To 
describe such a process we consider a combined system 
of the soliton and phonons 



^eff 



-ph 



^ph 



(16) 



where Lg = PX — H{P^ X) and the last term describes 
the low-energy phonons in the quadratic harmonic ap- 
proximation [28l [29], 

^2 h^n 



L 



ph 



dx 



(17) 



Here we use density-phase representation for slow bosonic 
fields il){x^t) = + t) exp(i(/?(x, t)) [30 . The cou- 
pling Lg-ph between soliton and low-energy phononic 
modes is given by the universal form, 

Ls-ph = -h^ ^{X, t) /tt - hN Lp{X, t), (18) 

where we have introduced displacement field ^{x^t) = 
TT p{y^ t)dy. The universal coupling ( 18) was derived in 
Ref. 2S| directly from the principle of translational and 
gauge invariance and was shown to remain valid away 
from weak coupling limit. It shows that it is the soliton's 
total phase and the number of particles ejected from the 
soliton core (or, rather their temporal change) that cou- 
ple to the phonon fields. The coupling (18) is also local. 



involving phonon fields at the position of the soliton. 
Dissipative dynamics of gray solitons. — The phonon 



action in Eq ( 17) is quadratic, thus we are able to directly 



integrate out the phononic degrees of freedom using the 
Keldysh formalism [26^ to obtain an effective action for 
the soliton. We obtain the following equations of motion 
modified by the phonons 

^-NdpN 



X = V 



2 



2/^ 



P = -dxH - 

he 



2 



h^ 



NdxN 



V 2c 



V 

2ne^ 



(19) 



(20) 



Here n = hn/me is a large number related to the Lut- 
tinger parameter K = ttk of the phononic Lagrangian 
(17) and V = dpH is the velocity of the soliton in the 



absence of phonons. These equations hold in the adia- 
batic limit, assuming slow deviation of soliton parame- 
ters. Similar expressions were obtained in Ref.[2T . 



As is well know from the theory of radiation in QED 



(see e.g. Refs.^[3T]) the dissipative Eqs. ([T9| and ([20]) 
are plagued with runaway solutions. The only consis- 
tent way to treat these equations is to solve for the soli- 
ton trajectory P(t), X(t), N{t)^ ^{t) in the absence of 
phonons and calculate the non-adiabatic corrections per- 
turbatively. As the number of particles N ejected from 



the soliton is conserved the r.h.s. of Eqs. (19) and (20) 



simplifies considerably. Turning to the energy dissipation 
rate, we have 



H = PdpH ^XdxH 



hhi 



2 c2 - V2 



(21) 



Using Eqs. (|6|,j7k one can show that ^ = -gNV /c^. 
Averaging Eq. ([2l]) over on e period using harmonic mo- 
tion V/c = vT^(7V72^cos(cjt/V^), we find the rate 
of slow change of energy of the soliton 



E 



1 

hnuo'^ 



Hdt 



1- ^- 



2hi<i 1 

2 



V^dt 



mg^N^ 
2^2 



N . 



(22) 



This allows us to calculate the life-time of the soliton. 



N 



mg 



N 



N'^ dN 



h^nuj^ Jo l-(7V/2/^)2 huu^ \2i^ 



where we have defined 



F{x) 



y'^ dy 



1. 1 

- log — 
2 ^ 1 



X . 



,(23) 



(24) 



For small x it behaves like F{x) while near x = 1 

it diverges logarithmically F{x) log — x). Using 

the fact that {N/2nf = E/E^^ we get the results stated 
in and after Eq. ([2|. 

Conclusions. — Hole-like Lieb II excitations can be gen- 
eralized for trapped systems by quantizing dynamics of 
gray soliton in non-uniform background. Despite their 
finite life-time, the obtained excitations are well defined 
for current experimental conditions. This can be verified 
by studying spectrum of phonon absorption, e.g using 
parametric driving proposed recently in Ref. [32^. To re- 
veal the quantum nature of the soliton its Heisenberg 
uncertainty AX = ^/h/SmuJN in the position has to 
be larger than its width, ^ = 2h^ /mgN [33 . Their ra- 

tio AX/^ = (0F/4)v/^I7?^ ~ {mg^/h^ujf'^Nll^ is 
already 0.14 for experiment in Ref. [18], and can be in- 
creased, e.g by reducing uo and making interactions g 
stronger by a tighter transverse confinement [5]. 
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